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I. INTRODUCTION 

The renormalization group (RG) is a set of tools and 
ideas used to investigate how the physics of an extended 
system changes with the scale of observation [iHS] . The 
RG plays a prominent role in the conceptual foundation 
of several areas of physics concerned with systems that 
consists of many interacting degrees of freedom, as is the 
case of quantum field theory, statistical mechanics and 
condensed matter theory (TJ |3HZ]- In addition it also 
provides the basis for important numerical approaches 
to study such systems f2| [HHTT] ■ 

Given a microscopic description of an extended system 
in terms of its basic degrees of freedom and their inter- 
actions, RG methods aim to obtain an effective theory, 
one that retains only some of these degrees of freedom 
but is nevertheless still able to reproduce its low energy 
(or long distance) physics. The effective theory is ob- 
tained through coarse-graining transformations that re- 
move those degrees of freedom deemed to be frozen at the 
observation scale of interest. For instance, given a Hamil- 
tonian H^^^ for an extended system one may aim to use 
successive RG transforms to obtain a sequence of coarse- 
grained Hamiltonians ,H^^\H^^\ . . ^ each an ef- 
fective theory describing the original system at succes- 
sively lower energy scales, c.f. Fig. [T] Broadly speaking, 
RG techniques fall into two categories depending on how 
the coarse-graining is implemented, namely momentum- 
space RG [6] and real-space RG [8HTT]. 

Momentum-space RG is applied to theories that are 
expressed in Fourier space. It works by integrating out 
high-momentum modes of a field and it is often asso- 
ciated to perturbative approaches. Instead, real-space 
RG is applied directly to theories that are written in 
terms of local degrees of freedom, say spins in the case 
of a spin system defined on a lattice. It is not linked to 
perturbation theory and can in particular be applied to 
strongly interacting systems. As proposed by Kadanov 
[1] , the coarse-graining transformation is implemented by 
replacing a block of spins with a single effective spin, a 
procedure refined by Wilson |5] and subsequently turned 



by White into the density matrix renormalization group 
(DMRG) algorithm [51 |S] , an impressively precise numer- 
ical tool to study one-dimensional systems. 

A major difficulty of momentum-space RG comes pre- 
cisely from the fact that it requires, as a starting point, 
a description of the system in Fourier space. Such de- 
scription is not always available and might not be ob- 
tained easily. Consider for instance a system of interact- 
ing spins, as specified by some generic spin-spin interac- 
tion. There, obtaining a momentum-space representation 
might be as difficult as solving the whole theory. In this 
and many other RG approach must be performed 

in real space. 




FIG. 1: Given a description of an extended system in terms of 
a Hamiltonian H^^'' one may use successive RG transforms T 

to obtain a sequence of Hamiltonians 

each an effective theory describing the original system at suc- 
cessively lower energy scales (or longer distances). Charac- 
terising the fixed points of the RG flow may help provide an 
understanding e.g. of the low-energy behavior of the system 
in the thermodynamic limit or of the stability of the system 
under various perturbations. 
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In spite of its indisputable success, the DMRG algo- 
rithm [HI [9] suffers from a shortcoming that has impor- 
tant implications. Because of the accumulation of short- 
ranged entanglement near the boundary of a spin block, 
the dimension of the Hilbert space used to effectively de- 
scribe the block must grow with each iteration of the RG 
transformation. As a result, for instance, unstable fixed 
points of the RG fiow (scale invariant critical systems) 
cannot be fixed points of the DMRG algorithm. Another, 
more practical consequence of this growth is that it lim- 
its the size of ID critical systems that can be analyzed 
and, most importantly, it severely limits the success of 
DMRG computations in higher spatial dimensions. 

Entanglement renormalization (ER) is a real-space RG 
method recently proposed in order to overcome the above 
difficulties [12]. The main feature of ER is the use of 
disentanglers. These are unitary transformations, lo- 
cally applied near the boundary of a spin block, that 
remove short-ranged entanglement before the system is 
coarse-grained. As a result, the effective dimension of 
the Hilbert space for a spin block can be kept constant 
under successive RG transformations, so that the ap- 
proach can be applied to arbitrarily large systems. The 
potential of ER — as well as that of the related varia- 
tional ansatz, the multi-scale entanglement renormaliza- 
tion ansatz (MERA) [T3] — to efficiently describe critical 
and non-critical ground states has been demonstrated for 
a number of spin and fermions models in one [12 [TH - 
[T9l \21\ and two spatial dimensions [22l - f29] . Most signifi- 
cant are the recent works in which ER has been applied 
to study a geometrically frustrated spin model [5S] and 
systems interacting fermions [2BVl29j : these are problems 
that cannot be addressed by quantum Monte Carlo tech- 
niques due to the sign problem. In addition, it has been 
shown that the MERA offers a natural representation for 
systems with topological order [30, 31 . Finally, several 
algorithms to compute the MERA have been put forward 

MM- 

In this work we explore the ability of entanglement 
renormalization to produce a sensible RG flow, one with 
the expected structure of fixed points and flow directions 
according to momentum-space RG, in Z? = 1,2 dimen- 
sional harmonic lattice systems. Such systems are an 
ideal testing-ground for ER. On the one hand, they have 
well studied properties [52H55] and can be fully character- 
ized in terms of correlation matrices, a fact that simplifies 
the analysis and conveniently reduces the computational 
complexity of ER calculations. On the other hand, an 
RG analysis of free-particle theories can be conducted 
simply and without approximations in momentum-space 
allowing for a comparison between the numerical results 
obtained using ER and the exact solution. The setting 
of harmonic lattices also allows for ER to be formulated 
in the language of bosonic modes, a formalism familiar 
to researchers in the areas of condensed matter physics 
and quantum field theory. 

The results of Sect. \V\ demonstrate that a real-space 
RG transform based upon entanglement renormalization 
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FIG. 2: An outline of this paper. The harmonic lattice Hamil- 
tonian Hrs, defined in terms of interaction between local de- 
grees of freedom, may be coarse-grained directly with a nu- 
meric implementation of a real-space RG transform, such as 



entanglement renormalization (ER), as outlined in Sect. IV 



Iterating the RG transform r times we get the r ^ effective 
Hamiltonian H^g- An effective Hamiltonian may also be ob- 
tained by first transforming the Hamiltonian to a momentum- 
space representation Hms , via Fourier transform of the canon- 
ical coordinates, as described in Sect. |ll] Momentum-space 
RG transforms may be applied (analytically) to Hms as de- 
scribed Sect. |III[ The dispersion relations from the real-space 
Hamiltonians H^g are compared to those of the correspond- 



ing momentum-space Hamiltonians, 

Sect, m 



H 



in the results of 



is able to reproduce the exact results from momentum- 
space RG to a high accuracy in D — 1,2 dimensional 
lattice systems, both for the critical and non-critical 
cases considered. Also demonstrated is the ability of 
the MERA to provide an efficient and accurate repre- 
sentation of the ground state of free boson systems, thus 
extending to the bosonic case the results of Ref. [22j . 
These results provide strong evidence that the ER ap- 
proach, which can be implemented without making use 
of the special properties of free-particle systems, could 
be used to investigate low-energy properties of strongly 
interacting systems not tractable with momentum-space 
RG approaches [TIHTB] . 

The paper is organized in sections as follows (see also 
Fig. [2]). Sect. [Tl| introduces the harmonic systems under 
consideration. In Sect. Ill the process of renormalizing 
the system in momentum-space is explained, highlight- 
ing conceptual features of the RG. The critical system 
and examples of relevant and irrelevant perturbations are 
considered. Sect. IV explains the details of the real-space 



RG implementation, both in terms of renormalizing the 
Hamiltonian and in terms of renormalizing the ground 
state directly. In Sect. |V]a comparison between results 
obtained from exact momentum-space RG and the re- 
sults from numerical real-space RG is presented. 



II. COUPLED HARMONIC OSCILLATORS 

In this work the low-energy subspaces of harmonic lat- 
tices in D = 1,2 spatial dimensions are to be analysed 
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FIG. 3: A comparison of an RG iteration for a D = 1 dimen- 
sional system (left) in terms of a dispersion relation, E(k), 
in momentum-space and (right) in terms of a lattice in real- 
space. The RG transformation maps the t^^ effective theory 
into a new effective theory H^'^+^'> while preserving the 
low-energy (or long distance) physics of the original theory. 
The three steps which comprise the iteration are described in 
detail in Sect. |III| for momentum-space RG and Sect. |IV A| 
for the real-space RG. There are many possible ways the real- 
space coarse-graining step (i) can be implemented, Fig. [4] 
describes two non-equivalent implementations. 



using both real-space and momentum-space RG transfor- 
mations. We begin with a brief introduction to harmonic 
lattices detailing equivalent representations in real-space 
and momentum-space coordinates, and the Fourier trans- 
form that shifts between such descriptions. For clarity 
the following derivations shall only be presented for the 
ID system as the generalization to 2D, or higher dimen- 
sional systems, is straight forward. The Hamiltonian for 
a chain of N harmonic oscillators each with mass m, an- 
gular frequency w, and coupled with nearest neighbors 
via 'spring constant' K, is written 



N 



-ff = Co > —p^ + -—Qr + K [qr+i - qr) 
^ — ' \ 2m 2 



N 



r=l 



(1) 



where in the second line we have chosen cq = 2m and 
defined K — mK for convenience. Note that periodic 



boundary conditions are assumed. The operators pi and 
Qi are the usual canonical coordinates with commutation 
[Pki Qi] — ifiSki- In our present considerations it is conve- 
nient to focus on the critical (massless) Hamiltonian 



N 

E 



(pl + 2K{qr 



(2) 



though the non-zero mass case will later be reintroduced 
in Sect. |III A| As a preliminary to the momentum-space 
RG, the Hamiltonian Hq shall be recast into momentum- 
space variables via Fourier transform of the canonical 
coordinates. The Fourier-space coordinates p and q are 
defined 
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r— 1 



1- = -7^2^'ir^ 

^ r—1 



-27rirK/N 



-—27rirK/N 



(3) 



Substitution of the Fourier-space coordinates brings the 
Hamiltonian of Eq. [2] into diagonal form, here a set of 
uncoupled oscillators 



{N-l)/2 



Ho 



E {pl 



K=-(Af-l)/2 



8i^sin^(^) ql 



(4) 



Defining k = 2TrK/{aN), with constant 'a' representative 
of the lattice spacing, the thermodynamic limit (N — oo) 
can be taken 



7r /a 



Ho 



p{kf + 8Ksm''{^)q{kf]dk. (5) 



The theory has a natural momentum cut-off A = ±n/a 
originating from the finite lattice spacing; taking the lat- 
tice spacing 'a' to zero recovers the continuum limit with 
corresponds to the field theory of a real, massless scalar 
field. Two equivalent representations of the harmonic 
chain have been obtained; that of Eq. [2] written in terms 
of spatial modes (amenable to numeric, real-space RG) 
and that of Eq. [5] written in terms of momentum modes 
(amenable to analytic, momentum-space RG). The dis- 
persion relation Eo(k) of the system, which describes the 
energy of momentum mode k and is known from the so- 
lution to a single oscillator, is given by 



Eo (k) = 2V2K\sm{ka/2) 



(6) 



The RG transformations of the Hamiltonian shall be cho- 
sen such that the resulting effective theory preserves the 
low-energy structure of the original theory, or equiva- 
lently, preserves the small k part of the dispersion rela- 
tion. As the low-energy part of the dispersion is gapless 
and linear we would expect that, under the RG flow, the 
renormalized Hamiltonians would tend to a fixed point 
that has a linear, gapless dispersion. 
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III. MOMENTUM-SPACE RG 

In this section the harmonic system of Eq. [5] which 
has been recast in a Fourier basis, is analysed using 
momentum-space RG. Although the RG analysis of in- 
teracting systems is more complicated that the free par- 
ticle analysis undertaken here, and often involves per- 
turbation theory, the basic procedure is the same. It 
is useful to consider the RG transformation as occur- 
ing in three steps, (i) Firstly the momentum cut-off 
is reduced, A i— A' = A/2, and modes greater than 
the cut-off are integrated out of the theory. As there 
is no interaction between the momentum modes of Eq. 
[sj this step is presently very simple; the cut-off is reduced 
A I— > A' = A/2 whilst leaving the form of the Hamiltonian 
for modes with momentum k < A' unchanged 



Tr/2a 



^0 = 



q{kY\dk. (7) 



A:— — -7r/2a 



(ii) Next the length associated to the system is changed, 
this can be implemented as a scaling of the lattice spac- 
ing^, a i-> a' = 2a, which gives 



TV I a' 



ka' 



q{kY]dk. (8) 



k— — -K/a' 



A change has been made to the observation scale of the 
system in terms of length. Next a change in the observa- 
tion scale is made in terms of energy. Indeed, in the final 
step (iii) the fields are rescaled 



1 



P (fc) p' {k) = -^p (k) 
q{k)^q'{k) = ^q{k).. 



(9) 



so that the new field operators have a modified commuta- 
tion relation [p'(fc), q'(fc)] = ih/2, in accordance with the 
desired change of energy scale. In principle the RG trans- 
formation is complete, however in this instance a further 
transform is required to recast the critical Hamiltonian 
into a manifestly invariant form. The field operators are 
rescaled once more 



P' ^^ p" = V2p' 



(10) 



^ The approach of rescaUng the lattice spacing is common in the 
context of condensed matter problems; equivalently we could 
have rescaled the momentum of the theory, k i—^ k' = 2k, as 
is the approach most often used for the RG in a quantum field 
theory setting. 



In contrast to the previous transform of Eq. [9j this trans- 
form is commutation preserving hence does not affect the 
physics of the system; namely the dispersion relation re- 
mains unchanged. Implementing the third step, together 
with the auxiliary transform of Eq. 
malized Hamiltonian Hq^^ is given 



10 the first renor- 



7r/a' 



p" (ky + 32K sin^ 



ka' 



h— — -K I a' 



q" {kf ] dk. 

(11) 

The RG transformation is summarized: (i) the degrees of 
freedom that are not relevant to the low energy physics 
are removed, followed by a changes of observation scale in 
terms of (ii) length [a ^-^ a' — 2a) and (iii) energy {h i— >■ 
h' = h/2). Starting from a theory with a natural length 
scale, the lattice spacing a, the RG transform is thus 
used to derive an effective theory with the new length 
scale a' . The scale factors chosen in steps (ii) and (iii) 
may depend on the implementation of the RG as well as 
the problem to which it is being applied. Iterating the 
RG transform (dropping the 'primes' from notation), the 
^th j-gnormalized Hamiltonian is 



■K j a 



p{kf + 2'^+^ K sin' [-^]q{kr ]dk, 



k— — TX I a 



with corresponding dispersion relation 



e'^J^ (k) ^ 2^+5 VK |sin (fca/2^+1) 



(12) 



(13) 



In the limit of infinitely many transforms, r cx), we 
get the Hamiltonian Hq°°'' at the fixed point of the RG 
flow 



n /a 



H, 



(oo) 



[p{kf + k{k^a^)q{kf^dk. (14) 



The fixed point Hamiltonian has a purely linear, gapless 
dispersion 



Ef'^ (k) = aV2k\k\ 



(15) 



as anticipated. 



A. Relevant perturbation 

In the previous section the massless harmonic lattice 
of Eq. [5 a critical system, was shown to be a (non- 
trivial) fixed point of the RG flow. We now consider 
the stability of this Hamiltonian under the addition of 
perturbations. Perturbations to the critical theory can 
be classified as being relevant or irrelevant depending on 
whether the deviations from the fixed point, induced by 
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the perturbations, grow or diminish under the RG flow. 
In order to study relevent perturbations a mass term Hj-d 
is reintroduced to the critical system of Eq. [5j as shall be 
shown shortly this term grows under the RG flow. Hence, 
it significantly modifies the low-energy physics from that 
of the unperturbed system. The mass term is diagonal 
in both real-space and momentum-space representations 



7T / a 



q (fc) dk. 



(16) 



k— — 7T/a 



The perturbed Hamiltonian H„i — Hq + m^H-cd is equal 
to the original Hamiltonian of Eq. [l] describing coupled 
harmonic oscillators of mass m. Since this perturbation 
term does not reorder mode energies {E{k) is still an in- 
creasing function of |fc|) the same RG transformations 
may be performed on H^c\ as on the unperturbed sys- 
tem Hq. This analysis shows that the r**^ renormalized 
perturbation term grows exponentially with the RG 
iteration r 



H, 



(r) 



rcl ; 



(17) 



Thus even the addition of even a small mass m to the 
critical system leads to a large difference between the 
perturbed Hm and original Hq Hamiltonians after only 
a few RG iterations; that is to say, the perturbed system 
tends to a different fixed point of the RG flow. This 
is further evidenced by consideration of the dispersion 

(t) 

relation Em of the perturbed system Hm 



(18) 



2^m 



O 



In the Fourier basis -ffirrci has representation 



7r/a 



Hr 



rcl 



sm' \ ] q{ky dk. (20) 



k — — n /a 



We consider the perturbed system Ha = Hq + aiJirroi 
with the perturbation strength chosen a > 0. The 
r"^ renormalized perturbation H^^J^^ may be obtained 
through the same sequence of RG transforms as applied 
to the unperturbed system 



■K j a 



H 



ir) 
rrcl 



22^+2 sin^ 



) q {kf dk 



k— — -K I a 



■K / a 



2-2-2 / [a^k' + 0{i^]]q{kfdK 



22r+2 



h— — -K I a 



(21) 



The perturbation Himi is thus exponentially suppressed 
under the RG flow. Equivalently, the addition of this 
term to the critical system has an effect on the low-energy 
physics that diminishes with each RG iteration, as is seen 
directly from the dispersion Ea of the perturbed system 

Hn 



4")(fc) -2^+1 



ka 



I2K + a sin"' 



ka 



= aV 2K \k\ + O (2-^^) , (22) 
which converges to the same linear, gapless fixed point 



as did the unperturbed system of Eq. 15 



(t) 

In contrast to the linear, gapless dispersion Eq of the 
massless theory described by Eq. [15] we now see a 
quadratic dependence of the energy on the momentum 
fc, together with an energy gap that grows exponentially 
with the number r of RG iterations. 



B. Irrelevant perturbation 

Perturbations that become smaller along the RG flow 
are termed irrelevant perturbations, as these do not af- 
fect the asymptotic, low-energy behavior of the system. 
In this section we construct and analyze an example of 
such a perturbation. The irrelevant perturbation Hind is 
constructed from neighboring and next-to-nearest neigh- 
boring quadratic couplings 



H 



N 

irrel = ^ ^ 



+ 1 



^ (9r+2 +qr f 



(19) 



IV. REAL-SPACE RG 

In this section an implementation of real-space RG 
based upon coarse-graining transformations of the lat- 
tice is described. Following the seminal works of Migdal, 
Kadanoff, and Wilson in real-space RG [TJ |2] , the coarse- 
graining transformation maps a block of sites from the 
original lattice C into a single site of a coarser lattice 
Let us consider a ID lattice C oi N sites, each site 
described by a vector space V. Wc divide C into blocks 
of two sites and, following Wilson, implement a coarse- 
graining transformation by means of an isomctry w 



(23) 



where is the vector space of two sites, V is the vec- 
tor space of a site in the coarser lattice C of A^' = N/2 
sites and ly is the identity in V. This coarse-graining 
transformation, which we shall refer to as a local projec- 
tion (LP) transformation, is depicted graphically in Fig. 
^ From an initial Hamiltonian H defined on lattice £ 
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FIG. 4: The coarse-graining step of the real-space RG, in 
which a block of two sites (each composed of M bosonic 
modes) of the original lattice C is mapped to a single site 
of the coarse-grained lattice C' , can be accomplished in dif- 
ferent ways. (Bottom) The simplest method involves applying 
a series of local projections, realized by isometric tensors w, 
whose function is to select and retain relevant block degrees 
of freedom. This shall subsequently be referred to as the local 
projection (LP) method. (Top) Entanglement Renormaliza- 
tion (ER) differs from the LP coarse-graining by the inclusion 
unitary disentanglers u, enacted across block boundaries, be- 
fore the truncation with isometries. 



we can obtain an effective Hamiltonian H' on lattice C 
via the transformation 



H' = W^HW, W = 



(24) 



The LP transformation has the property of preserv- 
ing locality of operators; for instance if the original 
Hamiltonian was a sum of two-body interactions, H = 
Sill then the effective Hamiltonian would remain 

a sum of two-body interactions, H' — X^i^i ^'a+i- 
alternative coarse-graining transformation, known as en- 
tanglement renormalization (ER), follows as a modifi- 
cation of the LP scheme. As with the LP scheme, we 
map two sites of C into a single effective site of C via 
an isometry w, however in ER one first enacts unitary 
disentanglers u 



r»2 



(25) 



across the boundaries of adjacent blocks, as show Fig. |4] 
Thus the effective Hamiltonian H' , as given by a trans- 
formation with entanglement renormalization, is 



H' = {W^'U'') H{UW), U = 



.®N/2 



(26) 



Inclusion of the disentanglers in the RG scheme, although 
having the unfortunate effect of increasing the computa- 
tional cost of numeric implementation, has profound im- 
plications regarding the ability of the method to produce 



a sensible real-space RG flow, as shall be demonstrated 
in the results of Sect. |Vl The action of the disentan- 
glers spreads local operators; if the original Hamiltonian 
was a sum of two-body interactions, H = X^iLi 
then the effective Hamiltonian would consist of three- 
body interactions, H' — X]i=i ''■i i+i j+2- However, un- 
der further transformations the interactions remain fixed 
at three body interactions, a property of the ER trans- 
formation referred to as a bounded causal width [13]. A 
more comprehensive introduction to real-space RG, in- 
cluding a discussion on the role of disentanglers and a 
description of ER schemes for 2D lattice systems can be 
found in Refs. (HHIIIS]. 

Starting from the orig inal Hamiltonian (/:(o),if(o)) ee 
(£, H) one could, by iteration of the RG map, generate a 
sequence of effective Hamiltonians defined on successively 
coarser lattices 

(27) 

with the transformation T representing either a LP or ER 



transformation. 



The LP transformation T^p 



is charac- 



terised by the corresponding isometry while the ER 
transformation is characterised jointly by an isom- 



etry and a disentangler, (w 



M „(r) 



). Equation 



27 



can be 



used to directly investigate how H changes under scale 
transformations allowing e.g. one to characterise the sta- 
bility of H under perturbations or investigate properties 
of the system in the thermodynamic limit directly. 

In order to make meaningful comparisons between ef- 
fective Hamiltonians defined at different length scales, 
say between H'^'^'^ and TJ^'^+i) defined on lattices £*-'^-' and 
£(•^+1) respectively, it is required that the dimension of 
the Hilbert space of a site in remains the sa,me as 

This 



23 



that of or equivalently that V = V in Eq. 
ensures that, for instance, two-body operators defined on 
and £("^+1) exist in the same parameter space, allow- 
ing a direct comparison between the coefficients that de- 
scribe the operators. Keeping the dimension of effective 
lattice sites in lattice L'^'^'^ constant between RG trans- 
forms is also desirable for computational purposes. In- 
deed, if the dimension of effective sites were to grow with 
each RG iteration then the computational cost of imple- 
menting the transformations would also grow, limiting 
the number of transforms which may be performed. An 
investigation of whether the LP and ER transformations 
can accurately coarse-grain Hamiltonians over repeated 
RG iterations while keeping a fixed local dimension is a 
focus of this work. 

Note that, in general, the disentanglers and isometries 
(u, w) that best preserve the low-energy space of a Hamil- 
tonian 

under coarse-graining will depend on the 
specific Hamiltonian H^'^'^ itself; in practice optimisation 
techniques [T4l - [l6] are required to compute the tensors 
(u, w). There are however some systems where analytic 
expressions for these tensors may be obtained [301 IS] • 
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A. RG of the Harmonic Lattice 

We now turn our attention to the analysis of the har- 
monic lattice system with ER^. Following the ideas dis- 
cussed in the previous section, several possible algorithms 
[HHTO] could be applied to study the low-energy sub- 
space of the harmonics systems directly. Application of 
these algorithms requires only that the Hamiltonian in 
question is composed of a sum of local interaction terms, 
as is the case with the systems we consider. However, 
as this study is focused on Hamiltonians containing only 
quadratic couplings, this property can be exploited in or- 
der to significantly reduce the cost of the numerical RG 
as well as simplify the analysis of the results. The Hamil- 
tonian of Eq. [l] describing a ID harmonic chain of N 
modes, can be written in a concise quadratic form 



H 



2N 

E 



(28) 



by defining a quadrature vector R — {p,q) with 



f Pi \ 

\PN J 



(29) 



where "H, henceforth referred to as the Hamiltonian ma- 
trix, is a 27V X 2N Hermitian matrix. The coarse-graining 
transformations shall be chosen such that the effective 
Hamiltonians also only contain quadratic couplings, de- 
scribed by some new Hamiltonian matrix "H'. Thus the 
RG analysis may be performed in the space of Hamilto- 
nian matrices H, as opposed to the (much larger) space 
of full-fledged Hamiltonians H and a more efficient real- 
ization of ER is possible. Retaining the quadratic form of 
the Hamiltonian entails limiting the disentanglers u and 
isometrics w which comprise the ER map to cannonical 
transformations, namely transformations that preserve 
commutation relations. Consider a transformation of the 
Hamiltonian matrix by a 2N x 2N matrix S 



n^n' ^ s^ns. 



(30) 



In order for the transformation to be commutation pre- 
serving it is required that the transform 5 be a sym- 
plectic matrix, S G Sp(27V, K). A symplectic transform 
can be characterized as leaving the symplectic matrix S 
invariant under conjugation, S'^'ES = E. Given our con- 
vention of grouping the quadrature vectors in Eq. 29 the 



^ The LP coarse-graining may be viewed as a simplification of 
ER in which the disentanglers u are set to identity transforms. 
Therefore only implementation of ER need be explicitly de- 
scribed. 



symplectic matrix takes the form 



(31) 



with Ijy as the N x N identity. Additionally, the Hamil- 
tonians under consideration take even simpler form than 
Eq. [28) as there is no coupling between p and q quadra- 
ture degrees of freedom in the harmonic chain, the Hamil- 
tonian may be expressed as 



H = fp 



q Hqq. 



(32) 



It is thus convenient to restrict symplectic transforms S 
to those which preserve the form of Eq. |32[ we only 
consider transforms of the type S — V (B V, with V a 
special orthogonal transformation, V S SO(A'^). It can be 
easily checked that ^ ® V is a symplectic transform. In 
fact, this is an element of the maximal compact subgroup 
of Sp{2N, M). The p-quadrature part of the Hamiltonian 
in Eq. [32] remains trivial under these transformations, 
allowing us to focus on the g-quadrature part Hq of the 
Hamiltonian matrix, which transforms as 



Ha 



V HqV. 



(33) 



Let us group a number M of contiguous bosonic modes of 
the ID harmonic chain together; each group of M modes 
shall henceforth be referred to as a site of the original 
lattice C. The disentanglers m, which act on two sites 
(that is, on 2M modes), are chosen as special orthogonal 
transforms u E S0(2M). Isometrics w are realized as a 
composition of a special orthogonal transform followed 



by a projection, w — wqW- 
Wo e S0(2Af), 



proj.7 



with 



(Oa/ e/M). 



(34) 



The transformation of the entire lattice is achieved by 
first constructing the direct sum of the local operators 



JV/2M 

W= . 



N/2M 



(35) 



Given the disentangler and isometry {u,w), the Hamil- 
tonian "Hq is coarse-grained into a new Hamiltonian H'^ 
defined as 



W'U' {Hq)UW. 



(36) 



By the definition of the isometry in Eq. [34] it is ensured 
that if the initial lattice C has M bosonic modes per 
lattice sites then the coarser lattices also have M 
modes per lattice site. As discussed earlier in Sect. IV 
keeping the number of degrees of freedom per lattice site 
constant between RG maps is necessary to allow mean- 
ingful comparison of operators at different length scales. 
Note that the number of modes per site M plays the role 
of a refinement parameter; choice of a larger M retains 
more parameters in the description of the effective the- 
ory, yielding more accurate results, at the cost of greater 
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entanglers and isometries as Eq. 36 Recall from Eq. 




2^' 



mum 



00 



{Hi) 



FIG. 5: An iteration of entanglement renormalization, broken 
into three steps, is depicted in terms of the direct sum struc- 
ture of the Hamiltonian matrix Hq. Dark shaded squares 
in Hq represent couplings within the site of AI modes, light 
shaded squares are the couplings between sites (at most next- 
to-nearest neighbor). Step(i), removing 'fast' modes is real- 
ized by transforming Hq by conjugation with the disentan- 
glers u and isometries w. Step(ii), rescaling the momen- 
tum, is achieved by removing the zero rows/columns from 
H'q. Step(iii), rescaling the fields, is achieved by directly scal- 
ing "Hq by a factor of 2. These three steps combined take the 
r**" renormalized Hamiltonian matrix Hq'^'' to the (r -|- I)*'' 
renormalized Hamiltonian matrix Hi^^^^ 



computational expense. The simplest choice of a one-to- 
one correspondence between bosonic modes and lattice 
sites, i.e. setting M — 1, does not give sufficiently accu- 
rate numerics, hence the need for grouping M > 1 modes 
into each lattice site. 

It is only for a proper choice of disentangler and isom- 
etry (u, w) that of Eq. 36 retains the low-energy 
subspace of the original Hq. Tms proper choice is found 
by optimisation over all possible (u, w). It is desired that 
the RG transform be optimized to project onto the mini- 
mum energy subspace of the original Hamiltonian; a ma- 
trix equation which describes the minimization can be 
written 



min(tr{ij;}) 



36 



(37) 



The ma- 



with the effective Hamiltonian H'^ as Eq. 
trix Hq describes a Hamiltonian that is translation in- 
variant between blocks of 2M modes, hence the trace of 
Hq (which describes the entire system) may be minimised 
by minimising the trace of an individual block of Hq. An 
optimisation method, based upon alternating updates for 
isometries w and disentanglers u, can be used to find suit- 
able (m, w) that minimise Eq. 37 and best preserve the 
low-energy space of the original Hamiltonian Hq. The 
optimisation method to be used here is similar to the 
general algorithm described in Sect. IV of Ref. fT6l. 

Assuming that suitable (u, w) have been obtained, the 
full real-space RG transformation of the Hamiltonian 
may be achieved in three steps, as illustrated Fig. [Sj 
analogous to the three steps of momentum-space pro- 
cedure described Sect. HI Firstly, (i) the Hamiltonian 



34 that an isometry w, which acts upon a block of 2M 
modes, consists of a special orthogonal transform fol- 
lowed by a projection, w = woWpmy The projection has 
form Wproj — (Om ffl Im) with the trivial (zero) part Om 
describing the M modes to be removed from the system 
and the identity part Im describing the M modes re- 
tained in the effective description. In the next step, (ii) 
the zero rows/columns, those which were acted upon in 
the previous step by Om, are removed from Hamiltonian 

Co)" 

matrix to form a new matrix Hq . (iii) The final step 
of rescaling the fields is realized by scaling the Hamil- 
tonian matrix by a factor, the same rescaling as Eqs. 
|9] and [TO] for momentum-space RG is realized by defin- 



ing H 



(1) 



2H 



^^'^ with "Hg^"* as the first renormalized 
Hamiltonian. Iterating this procedure r times gives the 
^th j-Qnormalized Hamiltonian Hq^\ In the thermody- 
namic limit, N ~> oo, the transformation can be iterated 
an arbitrary number of times to obtain a sequence of 
Hamiltonian matrices 



H ^ 

(38) 

each describing a theory with quadratic interactions be- 
tween at most next-nearest-neighbor sites, and each de- 
fined on an identical lattice C that has M bosonic modes 
per lattice site. 



B. Ground state RG 

In the previous section an implementation of real-space 
RG that could be used to coarse-grain harmonic lattice 
Hamiltonians was described. We now detail a similar 
procedure which allows ground-states of the harmonic 
lattices to be coarse-grained directly. Since we are deal- 
ing with systems of free-particles, the covariance matrix 
7 gives a complete description of the (Gaussian) ground 



state iV'Gs)- In the case of a Hamiltonian as in Eq. [32 
the convariance matrix is of the form 7 = 7p © 7^ , where 
7p and 7^ are defined as 



hq)ij = 2(V'Gsl9*gj IV'GS) 



(39) 



The derivation of analytic expressions for (7^, 7^) is a 
standard calculation and can be found e.g. Ref. |33) . 
Again, we choose the disentanglers u and isometries w to 
be canonical transformations. We obtain a sequence of 
increasingly coarse-grained states each defined by covari- 
ance matrix 7*^"^^ 



matrix Hq is transformed with direct sums of the dis- 



(40) 

with 7*^^^' = 7 as the original ground state. The coarse- 
graining transformations of the ground state shall be re- 
alized by symplectic transforms S acting in the space of 



9 



the covariance matrix, 7 7' = S'^jS. As there is no 
correlation between p and q quadrature coordinates, i.e. 
7 — Jp (9 jq, we may restrict consideration to symplec- 
tic transforms S S Sp(2A^, M) that are only of the form 
S = A® (A"^)'^ with A a real, invertible N x N invert- 
ible matrix. The covariance matrix 7 transforms under 
conjugation by matrix S, which implies that 



lq^iq = A-^(^q)(A-^f, 



(41) 



yeilding a new state 7' = 7p ® 7^. The disentanglers m, 
which act on two contiguous sites each of M modes, are 
thus realized as real, invertible 2M x 2M matrices that 
transform the covariance matrix as per Eq. |41[ Isome- 
tries w, which act on two contiguous sites within a block, 
are realised as w = woWproj with wq as a real, invertible 
2M X 2M matrix and "u^proj — (0^/ © Im) a projection 
onto M modes of the block. The direct-sum of the local 
operators (m, w) is constructed 



N/2L 

= 0^ 



N/2L 



N/2L 



(42) 



(43) 



in order to coarse-grain the entire lattice. Starting from 
the covariance matrix 7 = 7p © 75 describing the ground 
state of the original system and, for any choice of (u, w), 
a new state 7' = 7^ © 7^ can be obtained on a coarser 
lattice C through an ER transform 



I'p = {WproiW^U'^) Ip (C/M^oW^proj) , 

iq = {W^ro^W^'U-^) 1, {{U'TiW^TW^^o^) . (44) 

It is only for proper choice of tensors (u, w) that the above 
transformation correctly preserves the properties of the 
original state and produces a meaningful coarse-grained 
state. We now address the issue of how the proper ten- 
sors (m, w) can be found. In the application of the RG to 
the system Hamiltonian, both in momentum-space and 
real-space formulation, the modes truncated at each iter- 
ation were chosen as high-energy modes in order to leave 
the low-energy structure of the original theory intact. 
For the ground state RG a different criteria is required 
to judge which modes should be truncated from the sys- 
tem. The proper truncation criteria in order to preserve 
the ground state properties, proposed by White as part 
of his DMRG algorithm [8] , requires that the truncation 
of a block should be chosen to keep the support of the 
density matrix for the block. In the present formula- 
tion of bosonic modes, in which the representation of the 
state is given by a covariance matrix as opposed to a 
density matrix, this rule imposes that only modes in a 
block that can be identified as being in a product state 



with the rest of the system can be truncated and safely 
removed from the description of the state. Thus in com- 
parison with Hamiltonian RG, which was optimised to 
truncate modes such that effective theory had minimal 
energy, here we optimise for (u, w) so that the truncated 
modes have minimal entanglement with the rest of the 
system. 

The entanglement of a block with the rest of the system 
is known to be related to the symplectic eigenvalues of 7 
for the block [32l|33]. Let the covariance matrix (7)|b = 
{1p)\b © {'lq)\i3 describe the correlations within a block B 
of two A'/-mode sites. The 2M symplectic eigenvalues \i 
of the block B are the eigenvalues of the matrix formed 
by taking the product of matrices (7p)|b and (7q)|e 



\i = Spect{(7p)|e(7q)|e}- 



(45) 



The Heisenberg position-momentum uncertainty rela- 
tion, which here may be simplified as (p^) (9^) > 1/4, 
enforces that all symplectic eigenvalues are positive and 
have magnitude greater than unity, > 1 for all i. The 
entanglement entropy Si of a mode 'i' with symplectic 
eigenvalue is 



S, = 



f 



2 



/A[+l 
2 



(46) 



with f{x) = — xlogx. It is seen that the entropy Si is 
zero when mode i is in a minimum uncertainty state, 
Ai = 1, and an also that Si is increasing function of A;. 
It follows that, if a mode with eigenvalue Ai = 1 can be 
identified within a block B, then we are assured that the 
mode is in a product state with the rest of the system and 
may be safely truncated. Hence the tensors (w, w) should 
be optimized to minimize the eigenvalues Ai, and thus 
the entanglement entropy, of the modes to be projected 
out. The projection Wproj in Eq. 42 was defined to 
select the modes to be retained; we now construct the 
complimentary projector Wproj 



N/2M 



= (I2M - Wproj) 



(47) 



which projects onto the space of the modes to be trun- 
cated. The part 7 of the covariance matrix 7 that is to 
be projected out during the RG iteration may be written 
7 = 7p © 79 with 



Ip = {WproiWjU^^ 7p (uWoWp,, 

% = {WproiW^'U-') 7, {{U-T{WoTWp,.,i) • 

(48) 

Given that the modes projected out of block B are to 
have minimum entropy, the transforms (u, w) should be 
chosen to minimise the matrix equation 



min (tr{(7,)|B(7p)|B}). 



(49) 
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state into the tensors (u^*-*, ui^*-*) while preserving the 



TTN 




FIG. 6: (Top) Applying successive RG maps based upon ER, 
as depicted in Fig. |4j to the ground state leads to an ap- 
proximate representation of the ground state in terms of a 
set of disentanglers and isometries (u''^-', w''^') connected in 
a tensor network. This tensor network forms the multi-scale 
entanglement renormalization ansatz (MERA) [O], an ansatz 
for many-body states on a lattice, (bottom) Implementing the 
RG based upon an LP coarse-graining gives an approximation 
to the state in terms of a tree-tensor network (TTN) 36, 37 , 
a different class of ansatz for many-body states on a lattice. 



As with the energy minimization described by Eq. 37 for 
the Hamiltonian RG , this equation is optimised varia- 
tionally to find good disentanglers u and isometries w. 

The application of RG transformations to the Hamilto- 
nian could be interpreted as producing effective theories 
for the low-energy subspace of the original system. Is 
there a similar interpretation for the coarse-graining of 
the ground state? To address this question we assume 
that the ER map has been iterated r successive times on 
the original ground state 7^°^ to get the coarse-grained 
state 7*^"^^ . Each of the r coarse-grainings is characterised 
by a disentangler u and isometry w, thus the set of these 
transforms {u'^^\w''^^) for j = 1, 2, . . . , r characterise the 
sequence of RG maps. If the modes truncated during 
each iteration were in an exact product state, equiva- 
lently the eigenvalues of every mode 'i' truncated was 
Ai = 1, then the sequence of transformations could be 
inverted as follows. Starting from 7^^-', truncated modes 
(which were in a product state) are replaced back, and 
each of the transforms of Eq. |44] is inverted 



[u 



(--1) ,,,(--1) 



(- 

(50) 

as to recover the exact original state 7'"'. The coarse- 
graining of the ground state can be interpreted as stor- 
ing information about the short range properties of the 



long range information about the state; the set of ten- 
sors (it'^'^w^'^) together with state 7^'^^ thus serve as a 

representation of the original state 7^°^. The set 7'-'^^ 
and (u^*-*, w^*-*) form the multi-scale entanglement renor- 
malization ansatz (MERA) |13| , a variational ansatz for 
states on the lattice, c.f. Fig [6] [If the ground state RG 
is performed with an LP coarse-graining, as opposed to 
one based upon ER, a tree tensor network (TTN) |361l37j 
approximation to the state is obtained in terms of the 
isometries w^'^]. For the harmonic lattices we consider, 
as with most non-trivial models, the coarse-graining can- 
not be performed exactly and a MERA will be an ap- 
proximate, rather than exact, representation of the true 
ground state. In the present case, irrespective of how 
the transforms (m, w) are chosen, the modes to be trun- 
cated will still be slightly entangled as manifested in their 
symplectic eigenvalues, which will fulfill > 1. This 
entanglement will be ignored during the coarse-graining 
thus limiting the precision with which the original ground 
state 7^"-' may be recovered. 

In the present setting of free bosons, in which the ex- 
act ground state 7'-''' is already known, we compute the 
MERA and TTN approximations to the ground state via 
coarse-graining the ground state covariance matrix. 

The point of this exercise is that it allows us to assess 
the accuracy with which a MERA or a TTN can represent 
the many-body state. It also provides the opportunity to 
study the RG fiow of the ground state 7^"^^ under scale 
transformations. However, had our goal been to investi- 
gate properties of the unknown ground state of a system 
H with the help of real-space RG, it would have been 
absurd to assume knowledge of the exact ground state 
IV'Gs) from the beginning. In that case, an approxima- 



V'GS ) to the ground state may be found through 



tion 



e.g. variational minimization of energy ( V'gs 



H 



as 



as per Ref. [16; or through an alternative method P^llSj . 



V. RESULTS AND DISCUSSION 

A. RG flow of Hamiltonians 

In Sect. |III| the low energy subspace of harmonic chains 
were analysed exactly with momentum-space RG. This 
analysis yielded sequences of dispersion relations describ- 
ing the RG fiow towards a critical fixed point as well as 
the RG fiow resulting from adding relevant and irrelevant 
perturbation terms, as given by Eqs. 



15 18 and 22 



spectively. The same harmonic systems can also be anal- 
ysed numerically by applying successive real-space RG 
'transformations T, based upon ER, following the varia- 



tional method described in Sect. IV A This produces a 
sequence of Hamihonian matrices (TF'' ,n'-^\n'^'^\ . . .) , 
each an effective theory describing successively lower en- 
ergy subspace of the original system. By construction. 
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FIG. 7: Sequences of dispersion relations of Hamiltonians and 
after r = 0,1,3 RG transforms, comparing the real-space 
ER results (bold) with the exact momentum-space results 
(dashed, offset by 0.2). (Top Series) The critical Hamilto- 
nian Hq tends to a linear dispersion under the RG map with 
ER, in good agreement with momentum-space results. (Mid- 
dle Series) The addition of even small mass, m = 0.2, to 
the critical system gives a marked difference in the dispersion 
relation after r = 3 RG transforms. (Bottom Series) The 
effect of an irrelevant perturbation quickly diminishes under 
RG flow; by r = 3 iterations the original and the perturbed 
dispersions are virtually identical, Ea ~ -Eg . 
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FIG. 8: The entanglement entropy S of a site in the r**" 
coarse-grained ground state 7^'^' of the ID harmonic chain. 
(Left) In the critical (massless) regime the entropy of a site 
is infinite, as realized by the infinite constant il; however 
the change in entropy along the RG flow can be computed 
via a limiting process. The entropy of the state renormal- 
ized with an LP coarse-graining increases by a constant with 
each RG iteration, reproducing the logarithmic growth law, 
Sl = (1/3) logj L + c, as expected from conformal field theory 
[381 139) , whilst the entropy of the system renormalized with 
ER remains constant along the RG flow. Furthermore the se- 

.} rapidly 



quence of renormalized ground-states {7'^', 7'-^', 
converge to a fixed state 7* under the ER transforms. (Right) 
For several values of finite mass, the entropy of the states 
renormalized with the LP scheme saturate at a length scale 
governed by correlation length, whilst the theories renormal- 
ized with ER factorize into a product state (zero entropy) at 
the approximately same length. 



H^'^^ only contains quadratic couplings between nearest 
and next-to-nearest neighboring sites. The dispersion re- 
lations of the effective Hamiltonians V,'^'^^ are found by 
Fourier transform in a similar manner as presented for 
the original Hamiltonian in Sect. [H] Figure [T] shows the 
comparison of the analytic (momentum-space) and nu- 
meric (real-space) dispersion relations. 

The dispersion relations produced from numeric real- 
space RG with ER are seen to approximate the exact 
results to a high degree of accuracy for the three of the 
cases considered: a critical system, and relevant and irrel- 
evant perturbations on the critical system. Furthermore, 
the critical Hamiltonian converges to a fixed point 

of the RG flow, H^j'^ = H*, after r > 3 RG transforms. 
This fixed point Hamiltonian matrix "H* is manifestly in- 
variant under further transformations, T{'H*) = T-L* , to 
within small numerical errors. Consequently the isome- 
trics w and disent anglers u, which comprise the ER trans- 
form T, also converged to fixed points w'-'^' = w* and 



u^'^) = u* for r > 3. The Hamiltonian matrix T-La of the 
critical system with the addition of an irrelevant pertur- 
bation, Ha^ = Hq^'' +aH^^^J^-^, converges to the same fixed 

point as the unperturbed Hamiltonian "H^T'' — 'Hg^'' = H* 
for T > 3 transforms. 

B. RG flow of ground states 

On the other hand, in Sect. |IV B| we describe a real- 
space RG to coarse-grain the ground state of the ID har- 
monic chain. The coarse-graining transformations can be 
based upon either the LP or ER schemes of Fig. |4j The 
real-space RG produces a sequence of increasingly coarse- 
grained ground-states (7^°^ 7'^-'^^ 7^^^ . . .) , where 7'^"^ is 
the ground state covariance matrix of a harmonic chain 
with mass m. Fig. [8] displays the entanglement entropy 
S, as defined in Eq. |46j of a site in the coarse-grained 
state 7^'^'' as a function of the RG iteration r. 

The entropy of the ground state 7q^' of the critical 
chain, m = 0, when coarse-grained with the LP scheme. 
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FIG. 9: (Left) Sequences of dispersion relations (in non- 
rescaled energy units) for the (top) gapped and (bottom) crit- 
ical ID harmonic systems after r = 0,1,2,3 real-space RG 
transforms, comparing the results from the local projection 
(LP) method, -Elp', to those from the entanglement renor- 
malization (ER) coarse-graining, -E^p'. The numeric disper- 
sions produced by ER agree with the exact results obtained 
from momentum-space RG (dashed) as to be almost visually 
indistinguishable; though small errors are noticeable near the 
momentum cut-off, k = Ti/a. The dispersion relations given 
by coarse-graining performed with the LP scheme, whilst rea- 
sonably accurate after the 1^^ iteration, rapidly diverge from 
the exact result with successive RG iterations. In all cases the 
numeric RG was performed keeping Af = 4 modes per site. 
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Momentum Mode Number 

FIG. 10; Sequences of energy spectra (in non-rescaled en- 
ergy units) for the (top) gapped and (bottom) near-critical 
2D harmonic lattice system after r = 0,1,3 real-space RG 
transforms. The performance of the numeric real-space meth- 
ods, local projection (LP) and entanglement renormaliza- 
tion (ER), are benchmarked against the exact solutions from 
momentum-space RG (dashed). The energy spectra are ob- 
tained by sampling the dispersion relation E{ki,k2) on a finite 
grid of G points, with G chosen very large, and then ordering 
the values, {Ei < E2 < ■ ■ ■ < Eg}- The spectra of the sys- 
tems renormalized with entanglement renormalization, E'^^, 
retain a high level of accuracy through all r = 3 RG itera- 
tions, though do lose precision near the high energy cut-off 
of the effective theory. Energy spectra E^jJ of the effective 
theories obtained with an LP coarse-graining have diverged 
considerably from the exact result by r = 3 RG iterations. 
The numeric RG was performed keeping Af = 9 modes per 
site. 



increases by a constant with each RG iteration. More 
precisely, if we recall that a site in lattice £('^-* corre- 
sponds to a block of L = 2"^ sites of the original lattice, 
this growth of entropy reproduces the expected logarith- 
mic scaling, Sl = (1/3) logj i + c, for ID critical systems 
[38l [39] . This growth demonstrates that it is impossible 
for a critical ground state to be a fixed point of the LP 
scheme. 

The ground state of a system with finite mass m can 
also be analysed with the LP method. Initially, the 
ground state displays the same logarithmic growth of en- 
tropy as in the critical case, but it saturates approxi- 
mately after r = log2 C iterations of the RG transforma- 
tion, where C is the correlation length. 

Turning to the ER scheme, if disentanglers are included 
in the coarse-graining step then the entropy per site of 

(r) 

state 7q remains constant under RG transforms. This is 
made possible by the disentanglers, which remove short- 
range entanglement at each iteration. Moreover, for the 
critical case, the sequence of successively coarse-grained 
ground-states 7q^-' explicitly converges to a fixed point. 



7^'^) of a harmonic chain with finite mass m converged 
to a trivial fixed point, a product state 7**, under RG 
iteration. This occurs at the length scale of the correla- 
tion length C. The numeric results have been obtained 
by keeping M = 4 bosonic modes per lattice site. The 
MERA approximation to the ground state, obtained from 
successive coarse-grainings with ER, proves remarkably 
accurate. For the critical system, which is the hardest to 
analyse from a computational point of view, the MERA 
obtained from r = 9 RG iterations can reproduce the ex- 
act correlators of the ground state to error bounded by 
1 X 10-^. 



Local Projection vs Entanglement 
Renormalization 



7o 



7* for T > 3. The non-critical ground states 



Figs. |9] and 10 present numeric dispersion relations 
for the harmonic lattices in D = 1 and D = 2 spatial 
dimensions respectively, and compare results produced 
by coarse-graining with LP to those obtained by coarse- 
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FIG. 11: The mean energy -Brs, as defined Eq. 51 of the ef- 
fective theories obtained through real-space RG, based upon 
either an LP (sohd) or ER (dashed) coarse-graining, are com- 
pared against exact momentum-space results -Ems in terms of 
the relative error Ai5 — (-Ers — Ems) /-Ems- Performing the 
real-space coarse-graining with a larger number of modes per 
site M, which allows for more parameters to be kept in the 
description of the effective theories, is shown to give more ac- 
curate results for the real-space RG. (left) For the ID critical 
harmonic chain the mean of the energy spectrum obtained 
from renormalizing with ER, while keeping M — 4 modes 
per site, remains approximately 1% greater than exact value 
throughout the RG iterations. The spectra obtained from 
the LP method, even when using significantly larger M, only 
gave at best 10% accuracy after the same number of RG iter- 
ations, (right) For the 2D near-critical harmonic lattice, ER 
gives accuracies after 5 iterations of no less than 4% and 2% 
for Af = 4 and M = 9 respectively; this is in contrast to the 
LP method which is only accurate to within 60% for M — 4 
and 30% for M = 9. Barring small fluctuations, the accuracy 
of the coarse-graining with ER remains relatively stable with 
RG iteration. 



graining with ER. The dispersion relations for the ef- 
fective theories produced from the LP method are rea- 
sonably accurate after a small number of RG iterations. 
However, they rapidly diverge from the exact results 
in subsequent iterations. On the other hand, coarse- 
graining with ER is shown to keep significantly better 
support of the low energy subspace and, most impor- 
tantly, to maintain accuracy over repeated RG iterations. 
The numeric spectra obtained with ER for the 2D lat- 
tices of Fig. [To] displays a loss of accuracy towards the 
high-energy cut-off, which indicates difficulty in keeping 
a sharp cut-off numerically. However, this is of little con- 
cern as the primary interest lies in the low-energy physics 
and not the high-energy cut-off of the effective theory. 

The average mode energy E of an effective theory, de- 



fined in terms of its dispersion relation E(k), is 



(51) 



for a ID system and similar for 2D. The average mode 
energy is used to qualitatively analyze the accuracy of 
the effective theories obtained with real-space RG, as 



presented in Fig. 11 Keeping more modes-per-site M, 
hence more parameters in the description of the effective 
theory, increases the accuracy with which the numeric 
RG may be performed. However, even with the choice of 
a very large M, very large the LP method shows signifi- 
cant increase in error along the RG flow for both the ID 
critical and 2D near-critical harmonic systems. This fig- 
ure also confirms what was established visually in Figs. [9] 
and 10 that coarse-graining with ER not only produces a 
more accurate low-energy theory than with LP, but also 
maintains constant accuracy over successive RG maps. 



VI. CONCLUSIONS 

Real-space RG techniques, primarily in the form of the 
DMRG algorithm (SJ [3] , have proved an invaluable tool 
for the numeric analysis of low-energy properties of ex- 
tended quantum systems. However, this class of real- 
space RG methods (including the LP coarse-graining de- 
scribed here) suffer from an important deficiency: namely 
they do not reproduce proper RG flows. As demonstrated 
by Fig. |8j a ID critical system could not possibly be a 
fixed point of the LP map due to the growth of entropy 
along the RG flow, which is caused by the accumulation 
of short-range degrees of freedom. This deficiency in turn 
limits the size of ID critical and 2D lattices that can be 
analysed with such methods. 

In this work we have demonstrated, by comparison 
against exact results from momentum-space RG, that 
a coarse-graining based upon entanglement renormaliza- 
tion can reproduce proper RG flows. Demonstrations 
included the analysis of a critical Hamiltonian Hq and 
its ground state 70, which were shown to rapidly con- 
verge to non-trivial fixed points of the RG flow, H'q and 
7o respectively, and analysis of several non-critical sys- 
tems which converged to trivial fixed-points of the RG 
flow. By addition of an irrelevant perturbation to the 
critical Hamiltonian Hq, a new Hamiltonian Jia was ob- 
tained that described the same phase as Ho, but differed 
in the local interaction. As is required of a proper RG 
flow, "Hq, and Hq converged to the same fixed point Hq 
under the RG map defined by ER. It is important to 
note that the tensors {u,w)^ which comprised the real- 
space RG transform, were not chosen based on heuristic 
arguments or the desire for a particular outcome. In- 
stead they were found through optimisation based upon 



either energy minimization, as in Eq. 37 or attempting 
to retain the support of the ground state, as in Eq. [49] 

In addition, the coarse-graining transformation based 
upon ER was shown to induce a sustainable RG map; one 
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that could be applied arbitrarily many times without sig- 
nificant loss of accuracy and without the need to increase 
the local dimension M of the effective theories (so that 
the computational cost is also kept constant). The sus- 
tainability of the ER based RG map allows investigation 
of low-energy properties in arbitrarily large or infinite ID 
and 2D lattices, as has also been demonstrated in recent 
studies in which local observables are evaluated directly 
in the thermodynamic limit [1211211 US] and critical expo- 
nents computed [HI US] without the need for finite size 
scaling techniques. 

The implementation of entanglement renormalization 
presented in this work exploited properties of free- 
particle systems in order to reduce the computational 
cost. More general algorithms exist [HHTB] which allow 



implementation without making use of the special prop- 
erties of free particle systems. Thus it is possible to use 
real-space RG, based upon ER, as a means of investigat- 
ing the low-energy properties of strongly-correlated sys- 
tems where perturbative approaches are not valid. How- 
ever, since the implementation of ER in the interacting 
case is computationally expensive, especially so for 2D 
lattices, it remains unclear the level of accuracy that may 
be obtained in practice. Early investigations of 2D lattice 
systems with the full algorithm are promising p3ti28j . 
and the development and application of entanglement 
renormalization remains a subject of ongoing research. 
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